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ABSTRACT 

We investigate the possibility that the E p oc EJ relation between the peak energy E p of the vF v 
spectrum and energy output £ 7 for long-duration GRBs arises from the external shock produced 
by the interaction of a relativistic outflow with the ambient medium. To that aim, we take into 

account the dependence of all parameters which determine E p and £ 7 on the radial distribution of 

1/2 

the ambient medium density and find that the E p oc £J relation can be explained if the medium 
around GRBs has a universal radial stratification. For various combinations of GRB radiative process 
(synchrotron or inverse-Compton) and dissipation mechanism (reverse or forward shock), we find that 
the circumburst medium must have a particle density with a radial distribution different than the 
R~ 2 expected for constant mass-loss rate and terminal speed. 
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1 INTRODUCTION 

Lloyd, Petrosian & Mallozzi (2000) have established that 
the 25-1000 MeV fluence $ of bright BATSE Gamma-Ray 
Bursts (GRBs) is strongly correlated with the photon en- 
1 ergy E p ob3 ' > at which peaks the burst vF v spectral energy 
distribution. More recently, Sakamoto et al (2008) have 
shown that 83 Swift-BAT and HETE-2 bursts display a 
E p oha ^ oc ejj°- 52 ±° 11 correlation, with the burst fluence mea- 
sured at 15-150 keV, while Ghirlanda et al (2008) report 
E p obs) oc $°- 32 ± 05 for 76 bursts (with known redshifts), the 
burst fluence being calculated in the 1 keV-10 MeV range 
(i.e. bolometric). 

Lloyd et al (2000) found that the 8 GRBs with redshifts 
known at that time are not standard candles and, thus, the 
Ep" 1 " 1 - 1 — $ correlation is not due to cosmological effects but is, 
most likely, intrinsic. In that venue, Amati et al (2000) have 
shown that the intrinsic peak energy E p and the isotropic- 
equivalent burst output £ 7 at 1 — 10 4 keV are correlated, 
E p cc £°- 52 ±°- 06 j f or a se t of 9 bursts with known redshifts 
(most of which are among those used by Lloyd et al 2000). 
Later, Amati (2006) found that E p oc £°- 49 ± - 06 for a set 

of 41 GRBs, while Ghirlanda et al (2008) arrive at E p oc 
gO.54io.oi for 76 burstg ^ 

The lack of bursts with a high fluence and average/low 
peak energy bursts in the E p obs ' > — $ correlation is, evidently, 
not due to selection effects (i.e. at least half of that corre- 
lation is real), with the thresholds for burst triggering and 
measuring the peak energy possibly affecting only bursts 
with a high peak energy and low/average fluence. Ghirlanda 
et al (2008) and Nava et al (2008) investigate this possibil- 
ity and conclude that selection effects are negligible for pre- 



Swift bursts but do truncate the distribution of Swift bursts 
in the Ep^ — $ plane. However, as the range of peak ener- 
gies of Swift bursts is much narrower than that of the entire 
sample, they conclude that the E p oba ^ — $ is not an artifact 
of selection effects. 



2 POSSIBLE ORIGINS FOR THE 
E p oc £ 7 /2 RELATION 

1/2 

The simplest explanation of the E p oc £ 7 relation is that 
it arises from viewing geometry and/or relativistic effects. 
Such an explanation is generic, i.e. it does not make use of 
a certain mechanism for the production of the GRB. 

In the former framework, GRBs arise from narrow jets 
seen at various angles 9, the intrinsic burst emission being 
relativistically boosted by a factor D = [r(l — (icosO)]^ 1 ~ 
2/(r# 2 ), with P being the jet Lorentz factor and the view- 
ing angle 8 being larger than the both the jet opening 
and the relativistic beaming angle Y~ . Relativistic beam- 
ing of the comoving frame emission (denoted with primed 
quantities) implies that the observed burst peak energy is 
Ep — DE' p and the inferred isotropic-equivalent GRB out- 
put is £ 7 = D 3 £!y (the factor D 3 arising from D 2 for angular 

beaming and D for boost of photon energy). Hence, in this 

1/3 

scenario, the simplest expectation is that E p oc £ 7 , as- 
suming that the comoving-frame peak energy E' p and GRB 
output EL are universal (i.e. they have same values for all 
bursts) or, at least, uncorrelated. Toma, Yamazaki & Naka- 
mura (2005) obtain the E p oc £ 7 ^ 3 analytical expectation 
in a more sophisticated way (for a typical GRB spectrum) 
but their numerical integration of the Doppler-boosted emis- 
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sion yields E p oc £? t A for observer offsets that are compara- 
ble to (but larger than) the jet opening, which corresponds 
to higher energies E p and £ 7 . Using an annulus geometry 
for the GRB outflow (i.e. a hollow jet), Eichler & Levinson 
(2004) obtain a relation between the apparent E p and £ 7 

1 /2 

consistent with or shallower than the E v oc £ 7 relation. 

A potential problem with the off-aperture jet model is 
the expected distribution of GRB peak photon fluxes. The 
average photon flux (taken as a measure for the peak photon 
flux) of bursts seen at an offset angle less than 8 is larger 
than C(0) oc ^/(Ept-,) oc D 3 £! y /(DE' p t' J /D) oc D 3 oc 6~ 6 , 
where t 7 is the burst duration. The number of such bursts 
is N(< 8) oc 8 2 . Thus, the cumulative peak-flux distribu- 
tion expected in this model is N(> C) oc C _1/ ' 3 (for a 
volume- limited sample), which is flatter than that measured 
by BATSE (e.g. Pendleton et al 1996), showing N(> C) oc 
C _1 at peak fluxes between 1 and 10 photons/cm 2 s and 
N(> C) oc C~ 3 ^ 2 at peak fluxes above 10 photons/cm 2 s. 

Relativistic beaming of the emission from a jet wider 
than r~ and seen from a location within its aperture may 
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also be a possible origin of the E p oc £ 7 relation, as both 
quantities of interest, E p and £ 7 , are affected by the source 
relativistic motion. In this case, D ~ T, E p — TE' p , and 
£ 7 = r£ 7 (only one power of D because relativistic angu- 
lar beaming also reduces the source observed angular size 
by a factor D 2 relative to that of the entire source) , hence 
E p cc £j is expected if comoving-frame burst properties were 
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universal or uncorrelated. Thus, the E p oc £ 7 relation can- 
not be explained with just relativistic effects and requires a 
correlation of the comoving-frame peak energy E' p and GRB 
output £ 7 or a correlation of at least one of these quantities 
with the source Lorentz factor F (Schaefer 2003). Evidently, 
progress in this direction requires that a specific mechanism 
for the GRB emission generation is adopted (as done below) . 
In that venue, Zhang &: Meszaros (2002) showed that 
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the Ep oc £ 7 relation may be accommodated with internal 
shocks, by noting that the comoving-frame magnetic field of 
a Poynting outflow (or that generated through shock dissi- 
pation) satisfies B oc L^ 2 nt /RT, where L p / int is the out- 
flow's Poynting flux luminosity (or that of the dissipated, 
internal energy) and R is the radius at which the burst 
emission is produced. The GRB synchrotron emission peaks 
at E p oc 7 2 Br oc 7 2 L^ 2 nt /i?, where 7?n e c is the typical 
electron energy in the GRB source. Thus, one obtains the 
Ep cc £^/ 2 relations if (1) the outflow's Poynting (or internal 
energy) luminosity is a good measure of the GRB output (in 
the sense that the L p / int /£ 7 ratio is universal) and if (2) 7 
and R are universal (or not correlated with L p / int ). 

Note that the above argument applies to any dissipation 
mechanism. For internal shocks, the first condition above 
would lead to a constraint between the history of ejecta 
Lorentz factors (T(t)) and the distribution of ejecta mass 
with the Lorentz factor, while the second requirement for the 
GRB radius would constrain only T(t). As for the condition 
on the electron Lorentz factor, we note that, if electrons 
acquire a fraction e e of the outflow's internal energy, i.e. 
N^m e c 2 — e e U' , where N cc £^/Y is the electron number 
(£ k being the outflow isotropic-equivalent kinetic energy) 
and U' oc V'u oc V' L(RT)~ 2 is the internal energy (with 
V' — <hxR 2 ctL l being the volume of the GRB source), then 
7 oc e e L p / int / Lk, where Lt = £k/t-, is the outflow kinetic 



luminosity. The requirement that 7 is not correlated with 
Lp/int (leading to the E v oc £ 7 ^ 2 relation) implies that ei- 
ther e e oc L~j int or that Lu oc L p / int , otherwise one would 

obtain that E„ oc L°{ 2 . oc £~/ 2 . 

Similarly, constraints on some model properties are re- 
quired to explain the E p oc £^ 2 relation if the burst emis- 
sion results from Comptonization of the thermal radiation 
produced by magnetic reconnection or shock dissipation be- 
low the baryonic and/or pair photospheres. In this model 
(Meszaros & Rees 2000, Ryde 2004, Ramirez- Ruiz 2005), 
the GRB peak energy and luminosity are correlated be- 
cause both depend on the photospheric temperature. Rees 
& Meszaros (2005) have shown that, if dissipation occurs 

above the saturation radius, then E p cc T 2 L 7 1 ^ 4 , where L 7 
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is the GRB luminosity. Then the E p oc £ 7 relation re- 
quires a certain correlation of the burst luminosity with the 
photosphere's Lorentz factor. Within the same model for 

the burst emission, Thompson (2006) has shown that the 
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Ep oc £y' relation is obtained if the burst thermal radi- 
ation is produced at the stellar progenitor's photosphere, 
for which the rest-frame temperature of the photons is 
(L 7 /r 2 ) 1/4 , and assuming that the outflow opening 
is set by its lateral expansion at the sound speed (8 oc T _1 ) 
and that the collimation-corrected GRB output (oc L 7 # 2 ) is 
universal (as was first indicated by the analysis of Frail et 
al 2001 and later shown to be incorrect by that the GRB 
collimated output ranges over 2 decades - e.g. figure 1 of 
Ghirlanda, Ghisellini & Lazzati 2004). 

In this work, we present a possible origin of the E v oc 
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£ 7 relation related to the dynamics of the GRB source, 
assuming an observer located within the opening of the rela- 
tivistic outflow (i.e. viewing geometry is not at work). If the 
burst emission is synchrotron, then the peak of the GRB 
vF v spectrum is at E p oc -y 2 BT and the flux density at 
that photon energy is F p oc BFN, where 7 is the electron 
typical comoving-frame Lorentz factor and N the number 
of electrons in the GRB source. The GRB output being 
£ 7 ~ FpEptj, it follows that the E p oc £ 7 ^ 2 relation requires 

Br 7 2 oc BTj(Nt, f ) 1/2 . (1) 

If the burst emission were inverse-Compton scatterings of 
the synchrotron emission generated by same electrons, then 
Ep picks an extra-factor 7 2 and F p a factor r, the optical 
thickness to electron scattering of the GRB source. Then the 
Ep oc £ 7 relation requires 

BTj 4 oc Br 7 2 (iVt y T) 1/2 . (2) 

It is tempting to attribute the E p oc £ 7 ^ 2 relation to 
(i) variations from burst to burst of the BT factor, which 
appears both in the peak energy E p and the GRB output 
£ y , (ii) universality of 7, and (in) the remaining "dummy" 
parameters (N, i 7 , r) being either universal or uncorrelated 
with BT (so that they do not yield a different E p — £ 7 de- 
pendence). We note that variations in 7 (for synchrotron) 
or 7 2 (for inverse-Compton) from burst to burst that are 
larger than those of BT would induce a E p oc £ 7 correlation 
for either emission process. 

The burst duration t 7 , which is the only observable that 
appears in equations (1) and (2), has a spread of 1.5-2.0 or- 
ders of magnitude among long-bursts, which is comparable 
to the observed spread in GRB energy £ 1 at fixed peak en- 
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ergy E p (see figure 1 of Ghirlanda et al 2008) . This suggests 

1/2 

that the observed spread in the E p oc £ 7 ' relation requires 
some correlation among the dummy parameters, although 
it is possible that the range of measured £ 7 is smaller than 
the true spread because, for a fixed peak energy, bursts of a 
lower GRB output may fall below detection. 

If the GRB emitting electrons are accelerated at shocks, 
then it is unlikely that the product Br can vary among 
bursts while 7 is universal, because acceleration of electrons 
at relativistic shocks is expected to yield an electron Lorentz 
factor 7 that depends on that of the shock. As the latter 
bears a connection with the GRB source Lorentz factor V, 
a universal 7 requires either universal F and e e or an ad- 
hoc correlation of these parameters. In the former case, the 
1/2 

E p oc £y relation would rest entirely on variations in the 

magnetic field B among bursts. The nearly 3 decades spread 

1/2 

in observed E v and that E v oc B oc e s imply that the 
fraction eg of the internal energy stored in shock-generated 
magnetic fields has a range of 6 decades. Thus, a universal 7 
requires a mechanism for electron acceleration at shocks that 
is completely decoupled from to the generation of magnetic 
fields, which is an extreme requirement. For example, in the 
Weibel instability model of Medvedev (2006), proton current 
filaments created by the instability produce electric fields 

which accelerate electrons over distances of about the proton 
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plasma skin-depth, leading to e e ~ e B , hence the 3 decades 
range of observed peak energies E p would be associated with 
an electron 7 which is far from being universal. 

1/2 

Thus, it seems much more likely that E p oc £ 7 re- 
lation is not due just to variations in the BT term among 
bursts and that some or all of the other parameters appear- 
ing in equations (1) and (2) contribute as well. To include 
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their effect in driving the E v oc £ 7 relation, we assume 
that the outflow's energy is dissipated by shocks which ac- 
celerates electrons and generate magnetic fields that acquire 
quasi-universal fractions of the dissipated energy. Some jus- 
tification for the latter assumption is that, if the electron 
and magnetic parameters e e and e_g were correlated as for 
the Weibel instability model, then their variations among 
bursts would induce a E v oc £~/ 4 ' correlation for synchrotron 
emission and E p oc £ 7 ^ 6 correlation for inverse- Compton. 

In the following section, we study the implications of 
equations (1) and (2), representing the E v oc £ 7 / ' 2 relation, 
in the framework of external shocks. We note that this model 
has the drawback that the efficiency of the GRB emission 
should be small (below 10 percent) for those bursts with a 
large number (hundreds) of pulses (Sari & Piran 1997). The 
same can be done for internal shocks which, as discussed 
above, will lead to constraints on the distribution of the 
ejecta Lorentz factor with mass and ejection time (or vari- 
ability timescale) . A low GRB efficiency is also expected for 
this model (e.g. Kumar 1999). 



3 EXTERNAL-SHOCK EMISSION 

The external shock driven by the interaction of the relativis- 
tic ejecta with the burst ambient medium offers two pos- 
sible GRB sources: the reverse shock, which energizes the 
ejecta, and the forward shock, which sweeps-up the circum- 
burst medium. Denoting by F' the Lorentz factor of either 
shock as measured in the frame of the incoming gas (the 



ejecta or the ambient medium) , the shock jump conditions 
lead to an internal energy density in the shocked gas that is 
u — (r' — l)n'm p c 2 , where n' = (4r' + 3)no is the comoving- 
frame particle density in the shocked fluid and n' that in the 
unshocked gas. Thus, for a relativistic shock T' ^> 1), the 
typical electron Lorentz factor is 7 oc T' and the magnetic 
field is B oc Fn 1 ^ 2 , where n is the ambient medium density 
and r the laboratory-frame Lorentz factor of the shocked 
gas (i.e. the GRB source) , the latter being valid also for the 
reverse shock because the contact discontinuity between the 
two shocked media is in hydrostatic equilibrium, (i.e. pres- 
sure and internal energy density is the same behind both 
shocks) . 

For a source moving at Lorentz factor F, the burst du- 
ration is t 7 = R-y/F 2 , where R 1 is the GRB source radius, 
which results from either the spread in the photon arrival 
time across the visible area of angular opening F~ or from 
the observer duration of the source travel-time up to radius 
R (provided that the source is decelerating or accelerating 
slower than F oc R 1 '" 2 ). Adding that the optical thickness to 
electron scattering is r oc N/R 2 , the E p oc £l/ 2 relation of 
equations (1) and (2) becomes 

r' 2 rV /2 oc F'F(nNR-,) 1/2 (3) 
for synchrotron emission and 

r' 4 rV /2 oc r' 3 rjv(n/^) 1/2 (4) 

for inverse-Compton. 

Below, we investigate the conditions required for the 
synchrotron and inverse-Compton emissions from the re- 
verse and forward shocks to accommodate the E p oc £l/ 2 re- 
lation with the following simplifications: 

(1) the burst emission is produced before the reverse shock 
crosses the ejecta shell, i.e. before the deceleration of the 
external shock starts. One reason for this restriction is that 
the shock-crossing radius offers a "milestone" in the evolu- 
tion of the external shock that could be the GRB radius 7? 7 , 
while no such reference point exist after deceleration sets in. 
A second reason is that it would be unnatural for a decel- 
erating external shock to radiate episodically, once during 
the burst, until 100 s, and then again staring after 1000 s, 
during the afterglow, as observed in the X-ray emission of a 
majority of Swift GRBs (O'Brien et al 2006, Willingale et 
al 2007), 

(2) the entire emitting fluid moves at the same Lorentz fac- 
tor F(Rj) and is filled with the same magnetic field B(R 1 ), 
with the values taken at the radius were the burst emis- 
sion is released. As E p , F P E P , and F are power-laws in 
the shock radius, for a radially extended burst emission, 
their burst-averaged (E p — J E v dF p / J dF p ) and burst- 
integrated (j FpE p dt oc £j with dt = dR/F 2 ) values have 
the same dependence on i? 7 as their values at i? 7 , 

(3) the distribution with energy of the shock-accelerated 
electrons is softer than dN/dj e oc 7<T 3 above the typical 
7, so that the peak of vF„ is for the typical 7 electrons and 
not at a higher random Lorentz factor determined by elec- 
tron cooling and/or acceleration, 

(4) the electrons with the typical Lorentz factor 7 do not 
cool significantly during the GRB emission. We note that 
only a small of the BATSE bursts (Preece et al 2000) have 
the F v oc v^ 1 ' 2 spectrum below the peak energy E p ex- 
pected in the opposite case. If the 7-electrons cool during 



4 A. Panaitescu 



the burst, the vF v spectrum still peaks at the synchrotron 
or inverse-Compton energy corresponding to 7 (i.e. the peak 
energy E p remains unchanged), but the flux at E p picks a 
multiplying factor 7 c /7 oc (F'FnR-y)^ 1 for synchrotron and 
a factor (j c /~f) 2 for inverse-Compton, owing to that most 
electrons are at the cooling Lorentz factor j c oc F /(B 2 R 1 ) 
for which the radiative cooling timescale is equal to the burst 
duration. 

Thus, the treatment provided below is not sufficiently 
comprehensive and serves only as an illustration of the con- 
ditions required for the external shock to account for the 
E p oc 8-/ 2 relation. 

3.1 Dense ejecta (semi-relativistic reverse shock) 

For the evolution of the ejecta-ambient medium interaction 
at times before the reverse shock crosses the ejecta shell (i.e. 
before the standard deceleration sets-in), the shock jump 
conditions can be used to derive a fourth-degree equation 
for the Lorentz factor F of the shocked fluid (which is the 
same for both the shocked ejecta and the swept-up ambient 
medium). As shown by Panaitescu & Kumar (2004), the 
solution of that equation is 



R + , a self-consistent solution (s = 3.5) is found, for which 



To 



l + 2r (n/n ej ) 



1/2 



-1/2 



(5) 



where To is the Lorentz factor of the unshocked ejecta and 
n' e j their comoving-frame density. 

In the n' e j 3> 4F n limit (thin and dense ejecta shell), 
equation (5) leads to F < r , independent of the n/n' e j ra- 
tio, and to a mildly relativistic reverse shock of constant 
Lorentz factor F'. If F and V' do not change with radius, 
then equations (3) and (4) imply that the E p oc £^' 2 rela- 
tion is induced by a certain correlation of the ejecta Lorentz 
factor To with the radius i? 7 where the GRB emission is 
released. We focus on the forward shock emission because 
the mildly relativistic reverse shock is unlikely to yield an 
emission spectrum peaking in the hard X-rays. For the for- 
ward shock, r' = T ~ r and N oc nR 3 = R 3 ~ 3 for an 
ambient medium density stratified as n oc R~ 3 with s < 3, 
and N ~ const if s > 3. 

For synchrotron emission and s < 3, equation (3) re- 
quires Tq oc R 2 ^ 3 ^ 2 . If To were universal, this leads to an 
inconsistent solution s = 4. Thus To should vary among 

bursts, in which case the requirement imposed by the E p oc 
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£ 7 relation is that the GRB emission is released at a radius 
that is correlated with the ejecta Lorentz factor. Further in- 
vestigation can be done if R 1 is determined in some way. 
The termination shock of the wind expelled by the GRB 
progenitor is the only milestone expected in the evolution 
of the forward shock, though it is not evident how it could 
set the GRB radius; even that were achieved, the location 
of the termination shock should not be related to the ejecta 
initial Lorentz factor. Instead, we speculate that the loca- 
tion where the forward-shock GRB emission is released is 
tied to the radius R + oc (ft/To) 1 ^ 3- ^ at which the reverse 
shock crosses the ejecta (£ k being the ejecta kinetic energy), 
and after which the blast-wave is decelerated. Then and the 
E p oc £^ 2 relation can be obtained if s = 10/3, which is 
inconsistent with the starting assumption s < 3. 

For synchrotron emission and s > 3, equation (3) re- 
quires F 2 , oc R\ . Relating 7? 7 to the shock-crossing radius 



Ep OC T^ 7 7/4 ~ T-3pV4 



For inverse-Compton emission and s < 3, equation (4) 
requires F 2 , oc R 2 5 ~ 3 . If Fo were universal, the E p oc £l/ 2 re- 
lation would be accounted for by an ambient medium with 
s — 2.5. In this case, i? 7 = R + leads to _R 7 oc £ k 
and E p oc R^ 11 ^ 8 oc £ k '^ 2 . If r is not universal, the 
Ep oc £l/ 2 relation is obtained for s — 11/4, leading to 



E p oc f%r: 



-11/8 



oc F^£ 



17 c—ll/2 



For s > 3, the E„ oc £}/ 2 re- 



lation requires r ( 2 , oc R^ 1 ^ 2 , which for i? 7 = R + leads to 
s = 2.5, i.e. an inconsistent solution. 

Therefore, for a thin ejecta shell, the E p oc £^ 2 rela- 
tion can be explained with synchrotron emission from the 
forward shock if GRBs are produced at the radius where 
the reverse shock crosses the ejecta shell and if the ambi- 
ent medium around bursts has a universal n oc R- 3 - 5 radial 
structure, bursts with higher peak energy E p and GRB out- 
put _E 7 resulting for lower ejecta Lorentz factors To or larger 
ejecta kinetic energies £ k . The E p oc £^ 2 relation can also 
be obtained with inverse-Compton emission if n oc R~ 2 5 
and To are universal, bursts with higher E p and E-y result- 
ing from a lower £ k , or if n oc R~ 2 75 for all bursts if To 
is not universal, a higher E p and E 1 being obtained for a 
higher r or lower £ k . 



3.2 Tenuous ejecta (relativistic reverse shock) 

In the n' e j -C 4r ( 2 'ra limit (thick and tenuous ejecta shell), 
equation (5) leads to F = (r /2) 1/2 (n' ej /n) 1/4 < F and to 
a relativistic reverse shock with V' ~ ro/(2r) S> 1. Con- 
sidering that the radial width of the ejecta shell increases 
linearly with its radius, the comoving-frame ejecta density 
is n' e j oc (£k/Fo) I R 3 ■ Then, for a ambient medium with ra- 
dial density profile n oc R 3 , we obtain that the Lorentz 
factor of the shocked gas evolves as F oc R < - 3 ~ 3 ^'' 4 '. There- 
fore, if s > 3, the shocked gas motion is accelerated by the 
ram pressure of the incoming ejecta, starting from a value 
well below To (and remaining below it at all times). For 
s < 3, the GRB source is decelerating (but this deceleration 
is substantially slower than that after the reverse shock has 
crossed the ejecta shell) . 

In the following investigation, we drop the dependence 
of two quantities of interest, E v and E 1 , on the ejecta 
Lorentz factor To, i.e. we assume it to be universal, and 

determine the stratification index s that accommodates the 
1/2 

E p oc £j relation. In this case, bursts have different peak 
energies and GRB outputs because their emission is pro- 
duced at different radii i? 7 . If i? 7 is identified with the shock 
having crossed the entire ejecta shell, then the GRB radius 
is set by the ejecta kinetic energy and the duration of the 
ejecta release, which is about the same as the observer frame 



,1/(4- 



s ' (Panaitescu & Kumar 



burst duration: R + oc (£ k t- 
2004). 

If r were not universal, the E p oc £ 7 ^ 2 relation could 
be explained if the GRB radius i? 7 and To satisfy a certain 
relation. Then, relating R-, with the shock crossing radius 
will lead to a certain correlation among Fo, £ k , and t 7 , an 
avenue which we will not explore any further. 

The continuous injection of relativistic electrons (in 
the downstream region) of a Lorentz factor 7 e oc F' which 
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changes with the outflow radius will lead to an electron pop- 
ulation at Rj that has a power-law distribution with energy, 
dN/d'je oc 7 ( T 9 . The effective index q can be calculated by 
first determining the medium structure parameter s(q) that 
accounts for the E p oc El/ 2 relation, then the dynamics of 
the shocked fluid T(R) (which sets 7 e ) and the derivative 
dN/ dR of the electrons number, from where dN/d^ e can be 
obtained and the loop is closed to find the exponent q * . 
Because 7 evolves with R, one must check that the assumed 
location of the peak of vF v is consistent with the evolution 
of j(R) and the inferred effective index q of the electron 
distribution. 

3.2.1 Forward shock 

For the forward shock, r' = V oc B 7 S ~ 3)/4 and N oc R 3 ' s if 
s < 3, while N ~ const if s > 3. For s > 3, most of electrons 
have been accelerated before the GRB radius i? 7 and we 
have to find a self-consistent solution considering that the 
peak energy E p is either at 7(-R 7 ) (which we will denote as 
71) or at some electron Lorentz factor 70 corresponding to 
when relativistic electrons were first produced. We will make 
the simplifying assumption that 70 is a universal quantity. 

For a dN/d^e oc -y^ q electron distribution, with most 
electrons being at 70, the flux at photon energy Eo oc 70 Br 
is Fq oc FNB, while the flux at energy E\ oc J 2 BY is F\ oc 
TNB (70/71 (1) If the peak of vF v is at E , then it 

can be shown that the E p oc E^/ 2 relation requires s = 5, 
for which -y e oc R 1/2 , dN/dR oc R' 3 , hence dN/d-y e oc -y~ 7 . 
Given that d~f e /dR > 0, q — 7 implies that the peak of vF v 
is, indeed, at Eo, consistent with the starting assumption. 
For this case, E p oc i? 7 3 ^ 2 . 

(2) If the peak of vF v is at E\, then the E p oc £^ 2 relation 
leads to s = 3 + 4/(g + 3), 7 e oc R 1/( - q+3 \ q = 5 which, 
together with d^ e /dR > 0, implies that the peak of vF v is, 

in fact, at Eo, in contradiction with the starting assumption. 

1/2 

For synchrotron emission and s < 3, the E p oc £ 7 re- 
lation given in equation (3) is satisfied if s — 3.5, hence this 
is not a self-consistent solution. 

For inverse-Compton emission and s < 3, the E p oc 
El/ 2 relation of equation (4) requires s = 19/7, leading 
to 7 e oc iT 1/14 , dN/dR oc RT 5/1 \ and q = 5 which, to- 
gether with ofye/dR < 0, implies that the peak of vF v is, 
indeed, determined by the 7(-R 7 ) electrons. For this case, 

c, d-25/16 

E p oc ity . 

For inverse-Compton emission and s > 3, the flux at 
photon energy Eo oc 7 4 Br is Fo oc tFNB, while the flux at 
energy E x oc ^BY is Fi oc rTN B^o/^i) 2 ^^ . 

(1) If the peak of vF v is at Eo, the E p oc £^ 2 relation 
requires s = 1, incompatible with the assumed s > 3. 

(2) If the peak of f F v is at then the E p oc £ 7 ' /2 relation 
leads to s = 3 - 2/(g + 2), 7e oc ir 1/(29+4) , dJV/di? oc 
R~ q ^ q+2 \ from where q = 5 and s = 19/7, inconsistent 
with the s > 3 initial assumption. 

Therefore, the S p oc £^ 2 relation can be accom- 
modated with the synchrotron emission from the pre- 

* If, at each instant, a power-law distribution with index p is in- 
jected downstream, the cumulative electron distribution will have 
the index q if q < p and the index p if q > p. For simplicity, we 
assume that the former is always the case. 



deceleration forward shock if all GRBs occur in a n oc R~ 5 
medium, but at different radii, or by with the inverse- 
Compton forward shock emission if the ambient medium has 
a universal n oc R~ 19 ^ 7 stratification. In either case, bursts 
of higher E p and £ 7 are those occurring at smaller radii. 

3.2.2 Reverse shock 

For the reverse shock, V' oc oc _R 7 3 "'^ 4 and the 

number of energized ejecta electrons evolves as dN/dR oc 
R 2 (Ton' ej )(Po - (3), where fa and (3 are the lab-frame veloc- 
ities of the unshocked and shocked ejecta, respectively. For 
To > r, we have (3 - ~ (2r 2 ) _1 . Using n' ej oc iT 3 and 
r oc Tiy ^ , 0116 s,rriv6S at dN/dR oc i? (1 ~ s)/2 , from where 
N oc i? 7 3 ~ s '/ 2 for s < 3 and iV ~ const for s > 3. 
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For synchrotron emission and s < 3, the E p oc £ 7 rela- 
tion leads to a contradicting s = 5. For s > 3, assuming that 
the uF v spectrum peaks at Eo oc 70-Br, the i5 p oc El/ 2 re- 
lation requires s = 5, implying j e oc R" 1 ^ 2 , dN/dR oc i? 7 2 , 
thus q = — 1 which, together with d'je/dR < 0, implies 
that the peak of i/F^, is, indeed, at Eo- For this case, we 
obtain E p oc R^ 3 ^ 2 . For s > 3, assuming that the fi*^ 
spectrum peaks at E\ oc 7? Br, the E p cc £ 7 ^ 2 relation 
requires s = 3 — 4/(^ + 1), from where j e oc 
dN/dR oc i?( 3 -'?)/('?- 1 ), leading to q = 3 and s = 1, in- 
consistent with the starting choice s > 3. 

For inverse-Compton emission and s < 3, the B p oc 
E^ 2 relation requires 8 = 1, yielding j e oc R 1 ^ 2 , dN/dR = 
const, hence q — —1, thus d^e/dR > implies that vF v 
peak energy is determined by the 71 electrons. In this case, 
E p oc R-f . For s > 3, assuming that the vF v spectrum 
peaks at Eo oc 7oBr, the E p oc E}/ 2 relation requires s = 1, 
which is incompatible with the working condition s > 3. 
For s > 3, assuming that the vF v spectrum peaks at E\ oc 
7i-BF, the E p oc E~/ 2 relation leads to s — 3 + 2/g, implying 
7e cc ir 1/(29) and diV/di? oc R- 1 ' 1 ^, from where g = -1 
and s = 1, again incompatible with the starting condition. 

Thus, synchrotron emission from the reverse shock can 
account for the E p oc £ 7 ^ 2 relation provided that bursts 
occur at various radii in a n oc R~ 5 medium, while inverse- 
Compton emission can explain the same relation if n oc R^ 1 . 
For either radiation process, bursts of higher E p and E 1 are 
those occurring at smaller radii. 



4 CONCLUSIONS 

We have investigated the ability of the external shock (pro- 
duced by the interaction of relativistic ejecta with the am- 
bient medium) to accommodate the E p oc £ 7 ^ 2 relation be- 
tween the burst peak energy and isotropic-equivalent en- 
ergy release. First, we noted that it seems unlikely that the 
1/2 

E p oc £ 7 relation is due to variations of the quantity BT 
among bursts, with the electron Lorentz factor 7 being uni- 
versal and all other quantities (N, i 7 , r) in the right-hand 
sides of equations (1) and (2) not being correlated with E p . 
That is so because electron acceleration at relativistic shocks 
is likely related to the generation of magnetic fields and with 
the strength of the shock and because the number of radiat- 
ing electrons N and the burst duration t 7 could be related 
with the dynamics of the external shock. 

For that reason, we have identified the conditions that 
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lead to the E p oc £/ relation by taking into account all the 
quantities that determine E v and £ 7 . After making some 
simplifications (uniform magnetic field, single Lorentz fac- 
tor in the shocked fluid, negligible electron cooling), and 
considering only the external shock emission before the re- 
verse shock crosses the ejecta shell (as an interrupted burst- 
afterglow emission from the same decelerating outflow seems 
too contrived), we have determined the dependence of E p 
and £~t on the radius where (or up to which) the burst 
emission is produced, with allowance for both the reverse 
and forward shock, synchrotron and inverse-Compton emis- 
sions, and relativistic or semi-relativistic reverse shock. In 
the latter case, only the forward shock is expected to pro- 
duce the high-energy prompt emission, the E p oc £}/ 2 rela- 
tion requiring a correlation between the ejecta Lorentz factor 
and the GRB radius which does not have a plausible justi- 
fication. For that reason, we favour an explanation of the 
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E p oc £ 7 relation where the reverse shock is relativistic. 

In our treatment of that case, the burst emission is as- 
sumed to arise over a small range of source radii or up to a 
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certain radius, the E p oc £ 1 relation resulting from varia- 
tions in that radius from burst to burst. The reverse shock 
crossing the ejecta or the external shock encountering the 
termination shock of the progenitor's freely expanding wind 
are the milestones in the dynamical evolution of the reverse 
and forward shock, respectively, that could set the location 
where the burst radius is produced. This implies that the 
variations from burst to burst in the radius at which the 
prompt emission is released is due to either (1) variations 
among bursts in the kinetic energy of the ejecta or in the 
duration of ejecta release (for a reverse shock origin of the 
GRB), or (2) to the history of the mass-loss of the GRB pro- 
genitor shortly before its core collapse (for GRBs produced 
by the forward shock). 

Within the external-shock model for GRBs, the E v oc 
£}/ 2 relation can be accounted for by just the power-law 
radial stratification of the burst ambient medium density. 
For the four possible combinations of dissipation shock and 
radiation process, we find the following density profiles: n oc 
R^ 1 (inverse-Compton from reverse shock), 7?~ 19 / 7 (inverse- 
Compton from forward shock), and R~ 5 (synchrotron from 
either shock) . In general, a steep ambient profile is required 

1/2 

to explain the slope of the E p oc £ 7 relation because of the 
weak dependence of the source Lorentz factor on the density 
of the ambient medium (F oc n~ 1//4 - equation 5). 

None of the ambient medium stratification required by 
the E p oc £}/ 2 relation is the n oc RT 2 profile expected for a 
massive stellar GRB progenitor expelling a constant speed 
wind a steady mass-loss rate. Considering that the burst 
emission occurs at 10 16,J cm and that the wind termi- 
nation shock moves at ~ 10 km s _1 ), this implies that, in 
the last 1000 years before core collapse, the Wolf-Rayet 
progenitor of long-bursts had a varying mass-loss rate or 
wind terminal velocity. However, we do rule out that, by re- 
laxing the simplifying assumptions made here, the ambient 
medium density profile required to explain the E v oc £^ 2 re- 
lation with the external-shock emission becomes consistent 
with n oc RT 2 . 

As the burst model employed here is that of the external 
shock before the reverse shock crosses the ejecta (i.e. before 
deceleration begins) , the ensuing afterglow emission could be 



attributed to the emission from the reverse or forward shocks 
after deceleration, with allowance for injection of ejecta and 
energy after the burst, to account for the extended afterglow 
emission (if it is from the reverse shock) and the X-ray light- 
curve plateaus (if it is from the forward shock). Then, the 
general lack of continuity of burst-to-afterglow emissions, 
shown by the steep fall-off of the X-ray flux by 2-3 dex at 
the end of the burst, would lead to a rather contrived model, 
where the discontinuous burst-to-afterglow emission requires 
a temporary switch-off of the external-shock emission, fol- 
lowed by a much softer emission (the afterglows) . A simpler 
is that where the two emission phases, prompt and delayed, 
are attributed to different outflows, with the burst arising 
from a narrower jet whose bright, high-energy emission is 
produced only before the reverse shock crosses the ejecta or 
the external shock reaches the wind termination shock, but 
having a sufficiently low, collimated kinetic energy, so that 
its post-burst (forward shock) emission is dimmer than that 
from a wider, more energetic outflow producing the after- 
glow emission. 
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